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1. INTRODUCTION 


Recently proposed space structures have grown increasingly large 
and complex. These may be delivered to orbit by the space shuttle and 
then deployed/assembled on orbit. To reduce weight, efficient designs of 
such systems tend to lead to flexible, low-frequency, and often joint- 
dominated structures. Interaction between rigid body motion and 
structural deformation will likely occur. For efficient operation of 
system requiring component articulation, it is desirable to maneuver 
components as rapidly as possible. Operational speed is limited by 
excessive dynamic deformation if vibrations are not suppressed. In 
order to suppress excessive vibration response, active controls may be 
utilized. The control design is usually based on linear methods, however 
the articulation is governed by nonlinear equations, moreover, design 
methods use reduced structural models. To access these design 
performance as well as stabilities, analytical simulations are usually 
performed. 

Simulation codes for multibody systems such as DADS[1], DISCOS [2], 
and TREETOPS [ 3 ] use assumed mode approach to describe the structural 
deformations of components. This approach requires users to pick a 
reference frame for the component, discretize the component into finite 
elements, select component boundary conditions upon which the modes will 
be generated, solve the eigenvalue problem for deformation modes and 
select a modal set for the application at hand. Modal selection is often 
the most crucial part in the procedure. Since deformations of the 
component are defined as linear combination of the selected modes, the 
component can only deform in the space spanned by the selected modes . 
Therefore, results of the modal approach will be misleading if any modes 
are significantly excited were not selected. To predict which modes will 
get excited can be a difficult challenge in a flexible multibody system, 
since the system configurations are changing with time. Especially is 
thus true for many proposed future spacecrafts which have complicated 
geometry and joint hinges. 

To provide an alternative approach which circumvented some of 
these difficulties, the LATDYN computer code was developed for research 
purposes. The LATDYN program is finite-element -based. The user model the 
component with finite elements, instead of using truncated modes which 
have to be generated outside the multibody analysis codes. In order to 
separate the rigid body motion and small deformations in the finite 
element approach, a coordinate system is chosen to represent the large 
displacement and rotation of the element. Deformations of the element 
are then defined with respect to the rigid body configuration of the 
element. At the element-level, mass matrices are calculated. The 
component mass matrix is obtained by assembling each elemental mass 
matrices as is typically done in conventional small motion/deformation 
finite element methods. 

To form the system mass matrix, most multibody simulation codes 
impose nonlinear kinematic constraints on components that connect to the 
same joint. Instead of using constraint equations, the LATDYN program 
builds the ^linge degree-of-f reedom into the system equations of motion 
to connect components that share a common joint in a manner patterned 
after connectivity relations in conventional, small motion, finite 
elements. The mass of the interconnecting joint between the bodies 
represents a significant portion of the total mass and the orientation 


845 



of the joint's hinge lines play an important role in determining 
structural behavior. It is thus reasonable to construct the finite 
element program with the joints as a part of element connectivity. 


2. KINEMATICS OF BEAM ELEMENT 

The kinematics developed here is applicable to arbitrarily large 
displacement and rotational motion of a beam with small deformations. 
Consider the beam element with finite element nodes 1 and 2 at initial 
(undeformed) and current (deformed) configurations in an inertial X-Y-Z 
frame, as shown in Fig. 1. In order to specify the configuration of the 
beam element, it is necessary to define a set of generalized coordinates 
that uniquely define the global displacement of every point in the 
deformed element. For each node of the element, an xi-yi-zi (i-1,2) 
nodal reference frame having its x axis tangent to the neutral axis of 
the beam and y, z axes coincide with the principle axes of the beam 
cross section, is chosen to locate and orient the node in the inertia 
frame. Vectors ri (i-1,2) from the origin of the initial xi-yi-zi nodal 
reference frame to the origin of the current xi-yi-zi nodal reference 
frame define the global displacement of nodes 1 and 2. Transformation 
matrices Ti from nodal reference frames to the global frame define 
orientations of the nodal reference frame. Deformation of the beam and 
displacement of any point in the beam now can be determined using 

and Ti . 


Before deformations of the beam can be defined, rigid body motion 
of the beam has to be separated from the large displacement of the 
beam. It is chosen to specify rigid body motion of the element by use of 
a convected coordinate system X c -Y c -Z c whose origin is located at node 
1. Initially, orientations of the convected coordinate system coincides 
with the nodal reference frame of both nodes 1 and 2. As the element 
moves with large displacement and small deformation, the orientation of 
the convected coordinate frame, hence the rigid body motion of the beam, 
is determined by defining the Xc axis of the Xc”Y c -Z c frame always lie 
along the line connecting nodes 1 and 2, and the Y c axis to lie in the 
plane formed by the y axis of xi-yi-zi frame and the X c axis of X c -Y c -Z c 
frame. With these definitions, the convected coordinate system is 
uniquely determined. 

Deformations of the beam element are defined with respect to its 
rigid body configuration as 

Di = T c T Ti (1) 

where T c denotes the transformation matrix from the rigid body 
configuration (or the convected coordinate frame) to the global frame. 
Where Di is the difference in orientations between T c and Ti at any 
time step, due to flexural deformation, namely, a transformation from 
T c to Ti. Note that, Di can also be regarded as the transformation of 
Ti that rotates about a vector from the undeformed states to the 
current states. Assume that the rotation angles between Ti and Tq are 
small, then the components of this vector are the three rotation angles 
measured with respect to the three axes of T c [4] . Therefore, Di can be 
simply represented by 
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1 0zi $yi 

D H <(> 2 i i -<{> xi 

1_ ^ J 

where 4>yir and $ z i, are rotation angles of Ti about x, y, and z 

axis of T c , respectively. Physically, they correspond to flexural 
deformations of the beam element at nodes 1 and 2 . 

The rotation angles may be readily extracted from Di as follows: 

, T 


<|>xi = (Or 0, 1) Di (0, 1, 0) 

4>yi = (1, 0, 0)Di<0, o, 1) T 

<t» zi = (0, 1, 0)Di(l, 0, 0) T 

Substitution of Eq. 1 into Eqs.2-4 yields 

<t>xi = tc3 Tfc i2 
T 

<|>yi = tel fc i3 
4>zi “ t c 2 T til 


( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 
(7) 


where t c j and tij are 7 th column of T c and Ti, respectively. Since the 
x axis of T c lies on the line connecting nodes 1 and 2, the direction 
cosines of the vector from nodes 1 to node 2, which is the first column 
of T c , can be written as 


t c i = ( E 2 - ri + ro ) /I* 


(8) 


w h@re ro is a vector from node 1 to node 2 in the initial 
configuration. Since the Y c axis of T c lies on the same plane formed by 
tciand ti2 , the Z c axis of T c , hence the third column of T c , is the 
cross product of tci and tj.2/ namely. 


t c 3 = tclti2 


<9) 


The Yc axis can be easily obtained by taking the cross product of t c 3 
and t c i t 


t c 2 = t c 3 fc cl 

where a is the skew-symmetric matrix 


(10) 


0 -a ; 

a z 0 
a \j a y 


°y 

-ax 

0 
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with a x , ay, a z being the vector components of vector a. Note that 

a T = -a and that ab = -ba, which agrees with the vector product 
property axb ** -bxa . 

Let p = [x C/ y C/ z c ] be the position vector of an arbitrary point P 
in the beam element, defined with respect to T c . Define that 

£= x c /L 

T\ - Yc/ L 

C = Zc/ L 

where L is the length of the beam element. The displacement of a point 
P on a beam element due to flexural deformation may be expressed as 

u = NO (11) 

where N is a 3x6 matrix of shape functions similar to that used in the 
standard finite element method, namely 


N = 


-d-S> L C 

-(l-^)LTl 


(l-4^+3^ )LC (-1 + 4^- 3^ ) LTJ 0 

0 (£-2£ 2 +$ 3 )L -L^ 

(-^+2^ 2 -4 3 )L 0 


(-2^+3^)L(; (2^-3^ 2 ) : 

o (-^ 2 H 3 ) 
<^ 3 >l 0 


with T), and £ being x c /L, y C /L, and z c /L respectively, and <X>, the 
composite vector of rotation angles of nodes 1 and 2, 

T 

O = t ♦xl, ♦yl, $zl, 4^x2, 4>y2, $z2 1 

The total displacement of point P as shown in Fig. 1 is, in vector 

form. 


r p = r x + p + u - po 

where p is the vector [x Cf y Cf z c ] T and in algebraic form, 

r P = ri + T c p + T c u - T c0 p 0 (12) 

where T co is the initial transformation matrix of T c and P 0 is the 
initial position vector of point P in T c . Note that, axial deformation 
is implicitly included in the second term of the right-hand-side of Eq. 
12 . 

3. SUPER-BEAM ELEMENT 

In some multibody formulation, the joint connection between 
elements is imposed through constraint equations. Here, instead of 
introducing additional constraints, extra degree-of -freedom are added to 
the original element generalized coordinates to form a super beam- 
element consisting of joint and beam. 
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Consider a beam element with rigid joint bodies at both ends, as 
shown in Fig. 2. For simplicity, assume that no other element is 
attached to both joint bodies, and both joints have a one degree-of- 
freedom hinge. For each joint body, an Xi-Yi-Zi (i-1,2) body reference 
frame is chosen to locate and orient the joint body in the inertia 
frame. Vectors Ri (i-1,2) from the origin of the initial body reference 
frame to the origin of the current body reference frame define the 
global displacement of joint bodies 1 and 2. Transformation matrices Fj_ 
from joint body reference frames to the global frame define orientations 
of the joint body. Vectors that locate joint attachment points in joint 
bodies are denoted sj. (i— 1,2), defined with respect to joint body 
reference frames. Therefore, nodal displacements ri of the beam element 
can be represented by the joint body displacements Rj. as 


ri — Ri + Tisi Fio®i 


(13) 


To determine the relation between Ti and Ti, consider coordinate 

systems Xi'-Yi'-Zj. and xi -yi -zi that are located at a joint 
attachment point one fixed to the joint body and^one fixed to the nodal 

frame. Initially, let Xi -Yi -Zi and xi -yi -zi coincide with the z 
axes parallel to the hinge axis. Then, the difference in orientations of 
both systems at any time is the relative rotation ^ about the axis of ^ 


hinge. The transformation matrix from xi -yi -zi 
frame is 


frame to Xj 


0 [ = 


cos 0i -sin 0i 
sin 0i cos 0i 


0 

0 


0 0 1 


where 0i is the relative rotation angle. Denote Ti and Ti be the 

transformation matrices of Xi -Yi -Zi and xi -yi -zi frames with 
respect to Xi-Yi-Zi and xi-yi-zi frames, respectively. Then, the 
transformation matrix from joint body frame to the nodal reference frame 
can be obtained by sequential transformations as 


- Tiri ©i Ti 


(14) 


Substituting Eqs . 13 and 14 into Eq. 12, the displacement of an 
arbitrary point in the beam can be represented in terms of displacements 
and orientations of joint bodies at the ends of beam element, and 
relative rotations of joint degrees-of-f reedom, i .e . , 


r P = G(Ri, Ti, 0i, R2, ®2 ) 


(15) 


4. VARIATIONAL EQUATIONS OF MOTION OF A SUPER-BEAM ELEMENT 

The variational equations of motion of a beam element at time t f 
for a virtual displacement field that is consistent with the constraints 
is written as. 


849 


where 8r P is a virtual displacement of point P that is 

consistent with constraints, r p is the acceleration of point P, f p is 
body force density at point P , h p is the surface traction at point P, 

5e P is a kinematically compatible strain variation vector, T p is the 
associated stress vector at point P, and £2 and a are the volume and 
surface of the beam before it is deformed. 

By taking the variation of Eq. 15, the virtual displacement of 
point P is obtained as 

8r P = gR l 5Ri+gr l S^i + ge 1 S0i+gR 2 8R2 + gr 2 Sjt2+g0 J 502 {17) 

where 5rti is virtual rotation of H.. Therefore, the virtual displacement 
of a point P in the beam element is represented by the virtual 
displacements and rotations of the joint bodies, and virtual rotations 
of relative joint degrees-of-f reedom. 

The acceleration vector of a typical point can be obtained by 
taking two time derivatives of Eq. 15, which gives 

r p = g^Ri+grA+gejMgR^+grA+ge^ 

+gR 1 Hx + gr l Wi + ge 1 0i+gR 2 R2+gr 2 ®2+ge 2 02 <i8> 

where ©land 0) i are angular velocity and angular acceleration of the 
joint body i. 

Substituting Eqs . 17 and 18 into Eq. 16, the first term is 

- f ^8r pT r p dQ 


a 




Ri 

\ 



G>! 


= -[5Ri T , 8jt 1 T , 80! , 8 r 2 \ 8jc 2 t , 80 2 ]| 

M 

01 

r 2 

+ S (Ri,tO 1 ,0 1 ,R 2 ^2/02) ) 


i 

1 

o : e- 
1 

1 


(19) 

where M is the generalized mass matrix, 



M = I H 

Ja 


gRx T gRi gRi gn 
gr, T 9n 


gR^gex 

grx T 90x 

gex T gex 


gR, T gR2 

grx T 9R2 

go, T gR 2 

gR 2 T gR 2 


gRi T gr 2 

grx T ^r 2 


gox T 5r 2 

gR 2 T gr 2 


9r 2 T 9r 2 


Symmetric 

and S is quadratic in velocity. 


S = 


t 


gRx T (gR 1 Ri+gr l «i+g0 i ei+g R2 R2+gr 2 t02+g0 2 02) 
g ri T(g Rl R 1 +gr 1 O)i+g0x0i + gR 2 R2+gr 2 O)2+g0 2 02) 
g 0i T (gR 1 Ri+griC0i+g0,0i + gR 2 R2+gr 2 c»2+g0 2 02) 
gR 2 T (gR 1 Ri+griW 1 +g0 1 0i+gR 2 R2+gr 2 ®2+g0 2 02) 
gr 2 T (gRiRi + gri®i+g0i0i+gR 2 R2+gr 2 o)2+g0 2 02) 

L g© 2 T (gRiRi + griWi+g0x0i+gR 2 R2+gr 2 ©2+g0 2 02) 

Similarly, the second and third term in Eq. 16 become 


f 8r pT f p dQ + I 8r pT h p do 

J a Jo 


= [8R 1 T / 87i 1 T ,50 1/ 5R 2 T / 67C 2 T ,8e 2 ]Q 
where Q is the external generalized force vector. 


s-|J 


gRx TfP 


g R x T h p 

grx TfP 


grx ThP 

g 01 TfP 

dQ + f 

g 01 T h p 

g R 2 T f p 

Jo 


gr 2 TfP 


gr 2 TfaP 

_ 


g© 2 T h P _ 


da 


gRx T g© 2 

grx T 9o 2 


g0x T 90 2 

gR 2 T g© 2 

gr 2 T g© 2 

g0 2 T 9o 2 


dn 


( 20 ) 


For a Bernoulli beam, the right hand side of Eq. 16, or the 
virtual work done by the internal force, may be expressed as 


I 


8 e p Vd £2 


dQ 
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( 21 ) 


= I (EAuiSux+EIyU^Uj+EIzUjSuj'+GJ^x^xJdx, 

Jq 


= [ 6 R 1 T / 57 t 1 T / 50 lf 5 R 2 T , 57 c 2 T , 50 2 ]u 

where ui is the axial deformation, 

Ui«V (r 2 “r 1 +r 0 ) T (r 2 -r!+r 0 ) - Vr 0 T r 0 

U 2 and U 3 are bending displacements in the y and z direction of the 
neutral axis, which can be obtained from Eq. 11, and U is the 
generalized internal force vector. 

Substituting Eqs . 19-21 into Eq. 16, the variational equation of 
motion of a super-beam element can be written as 


[sRx , bn 1 r 50x, 5 r 2 , 57c 2 , 50 2 ]| 





Ri 



M 

0i 


r 2 


Ct> 2 


_ 02 - 


+ S(R lf CO 1 ,0 1/ R 2 ,CO 2 ,e2) + U - Q 


for all virtual displacements 5 Ri and 5R2, virtual rotations 5 tti and 
5 j 12 , and virtual hinge rotations 5©i and 502 that are consistent with 
the constraints. 

5. APPLICATIONS 

5.1.1 Analysis of A Deployable Space Structure 

The deployable space structure shown in Fig. 3 is a 20 meter long, 
triangular cross section, joint dominated truss structure, referred to 
as the Mini-Mast. The structure is used at NASA Langley as a ground 
test article for the development of research techniques in structural 
dynamic characterization of large space structures and control of 
flexible structures. A total of 18 bays, each 1.12 meter long, make up 
the 20 meter length of the beam above the deployer mechanism. Figure 4 
illustrates two deployed bays of the beam design in more detail. One 
bay of the truss beam consists of three longerons, three diagonal 
members and a batten triangular truss whose cross-section fits inside a 
1.4m diameter circle. The longerons and diagonal members are connected 
to batten triangles at each corner (three corner bodies are built into 
each corner of the batten triangle) by revolute joints. A sketch of a 
corner body is shown in Fig. 5, primarily to indicate geometric 
complexity . 

The system is deployed/retracted with two bays at a time. During 
deploying/retraction, the vertices of two batten triangles are held 
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fixed in orientation while a third batten triangle, located between the 
two fixed ones, rotates about the longitudinal axis. Upward/downward 
forces are then applied to deploy/retract the system. Revolute joints in 
corner bodies at each apex of the triangular cross-section of each bay 
and a nearly over-center hinge in each diagonal member allow the beam 
truss to deploy/fold into a repeatable beam/stack, as shown in Fig. 4. 

At the final stage of deploying, mid-hinge of the diagonal member is 
locked up to ensure that the system becomes a structure. It is reopened 
as retraction starts. Since each two deploying bays are symmetric to the 
middle batten triangle, it is sufficient to analyze only one bay of the 
system. 

The objective of the analysis is to determine loadings on flexible 
members during deployment of one bay of the truss beam. Because of the 
symmetric geometry of the system, corner bodies are constrained to move 
on a 1.4m diameter cylinder during deployment. The longerons and 
diagonal members are deformed, due to kinematic constraints imposed at 
the joints, during deployment. Therefore, orientations of revolute 
joints of the longeron and diagonal member play an important role in 
design of the truss beam. A set of properly designed revolutes will 
decrease deformations, hence decreasing the force required to deploy the 
truss beam. The system is designed such that it is not deformed in its 
fully retracted position and this then serves as a good starting point 
to analyze the system response during deployment. 

5.1.2 LATDYN Model Description 

The model can be simplified by taking the advantage of symmetric 
geometry of triangular cross section of the system. By constraining the 
upper triangle to only move along and rotate about the longitudinal axi3 
of the truss, the LATDYN model of one bay of Mini-Mast reduces to 3 
flexible longerons, 3 flexible diagonal members, and 2 batten triangles 
that are modelled as rigid bodies, as shown in Fig. 6. The lower batten 
triangle, batten triangle A in Fig. 6, is grounded. The upper batten 
triangle, batten triangle B in Fig. 6, is driven up and down. The 
batten triangles are connected to the longerons and diagonal members 
with revolute joints at each corner, respectively. The geometry and 
material properties of the longeron and diagonal members are listed in 
Table 1. Initial configuration of the model is chosen with the system 
in its fully packaged position, as shown in Fig. 6. Locations and 
orientations of each revolute joints of longeron and diagonal member are 
listed in Tables 2 and 3, respectively. 

5.1.3 Results 

The system is deployed by driving the upper triangle in the 
longitudial direction without constraining its rotation about the 
longitudial axis. The driving constraint is 

Z=L [t— 2-sin(^-) ] , t<T 

T 27t T 

Z=L, t£T 

where L is length of the longeron, T is total deploying time, and Z is 
the height of the upper triangle. The deployment moves the upper 
triangle a distance L in the z direction in T seconds. In the 
simulations which follow, T is taken as 1.0 second. 
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Figures 7-8 show the variations of the bending and twisting 
moments of the longeron at its midpoint, with the z displacement of the 
upper triangle. Figures 9-10 show the bending moments of diagonal member 
A and B while figures 11-13 display the time history of the axial forces 
of longeron and diagonal member, showing that all truss members in the 
system are actually in compression after deployed. This tends to 
increase bending stiffness of the truss beam. Figure 14 shows the LATDYN 
results for the bending moments of the longeron at the end joining the 
upper triangle. Also shown are the predictions of the Astro Co. [8]- the 
original Mini-Mast designer and producer, and the best result authors 
can achieve using the assumed mode approach. In Fig. 14, both LATDYN and 
Astro results predict a minimum bending moment in the longeron when it 
rotates about 50 degrees which is reasonable from the geometry of the 
structure. Results from the assumed mode approach indicate that the 
deformation modes used in the analysis are not complete, which is not 
apparent without another analysis test results. 


5.2.1 Flapping Motion of Rotor Blade 

A number of problems arise which make it necessary to study the 
effects of flexibility on blade motion. For example, the affect of 
flexible motion on the performance, stresses occur in the deformed 
blade, and interactions between the rotational speed and the natural 
frequencies of the flexible blade. An aditional complicating fator is 
that due to the stiffening effect of centrifugal force, natural 
frequencies of the blade increase as blade rotation speed increases. 

A simplified model of an articulated blade with no flap hinge 
offset or spring restraint is shown in Fig. 15. Initially, the blade is 
straight and tilted 0.157 radian. With no initial hinge velocity, the 
blade rotates at a constant speed. Since the centrifugal force always 
acts radially outward in a plane normal to the rotation axis, it acts as 
a spring force opposing the blade flap motion, hence initicating the 
flap motion and deformation. 

Simulation of the same flapping blade, using an assumed mode 
approach, producing a diverging solution as is reported in Reference 5. 
This is because the geometric stiffening effect is not properly 
accounted for [6] [7] . 

5.2.2 LATDYN Model and Results 

In the simulation of the flapping blade using LATDYN, Rotation 
speed of the blade is kept constant in each simulation and gradually 
increased in succeeding simulations, starting with 1 rad/sec and going 
up to 9 rad/sec. Frequencies are calculated from the transient response 
of the simulation using a Fast Fourier Tranform. 

Figures 16 and 17 show the bending moment of the blade at the 
midpoint when it rotates at 3 and 6 rad/sec. The results clearly 
indicate that the natural frequency of the first flapping bending mode 
increases as the rotation speed increases, due to the centrifugal 
stiffening effect. Figure 18 displays natural frequencies of the first 
bending modes for different rotation speeds, compared to the solutions 
derived by Southwell [9] . Good agreement between the LATDYN results and 
the Southwell solution is shown: Also shown (dotted lines) in Fig. 18 
are different harmonics of the rotor speed. As shown, the natural 
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frequency of the first mode intersects with the third harmonic around 8 
rad/sec, fourth harmonic around 4 rad/sec , fifth harmonic around 3 
rad/sec, and so on for higher harmonics. A resonance may then occur when 
the blade rotates around these speeds. Figure 19 shows the bending 
moment of the blade when it rotates at 8 rad/sec, which show that the 
magnitude of the bending blade keeps increasing with time. The frequency 
of the blade is about three times the rotational speed. The magnitude of 
the response in Fig. 19 may not increase indefinitely, but may represent 
a beating phenomenon with the period of the beat depending on the 
closeness of 8 rad/sec to the intersection point. 

CONCLUSIONS 

A three dimensional, finite element based simulation tool for 
flexible multibody systems is presented. Hinge degrees-of-f reedom is 
built into equations of motion to reduce geometric constraints. The 
approach avoids the difficulty in selecting deformation modes for 
flexible components by using assumed mode method. The tool is applied to 
simulate a practical space structure deployment problem. Results of 
examples demonstrate the capability of the code and approach. 
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Tali 1 e 1 Ccunii'liy iinJ H»l«*ilol l»t t»|ini ll ins 


I 

un i l> 1 ninct i*i 

Out <r f> i uiii"l i-i 


Longeron 

O.OH91 m 

0.02019 in 

hliiq. Member 

0.0U15 m 

0.01511 in 

M< win 1 nn of Elasticity 1.1 ‘IK* 11 ll/m 2 

Shear of Elasticity - 1 . 0E_» 9 II /in 2 

v * I.616EO kq/m^ J 

Table 2 Joints I.iiciit Ion 


Joint Name 

X 

Y 

Z 

Longeron A1 

-0 . 60598 

-0. 34987 

0.02379 

Longeron B1 

-0,23372 

0.65954 

0.02379 

Diagonal A1 

-0,58011 

-0.34987 

0.02379 

Diagonal B2 

-0.49966 

-0.50693 

0.02379 

[Manorial AIM/ 

-0,14861 

0.26657 

0.023/9 

Table 

3 Joints Orientations 

Joint Maine 

X-component 

Y- component 

Z-componont 

Longeron A1 

0.91838 

0.07838 

0.38786 

Longeron B1 

0.64757 

-0.65591 

-0.38786 

Diagonal A1 

0.98896 

0.09973 

-0.10956 

Diaqonal B2 

0.88133 

-0. 39542 

-0.2586B 

Diagonal A1B2 

-0.88076 

0.47339 

0.01163 




Figure 1 Kinematics of Beam Flement Figure 2 Kinematics of Super-Beam Element 
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Figure 3 * Deployable Space Structure 


Figure < 


Longeron goes Diagonal goes 

in Ibis slot in this slot 



Figure 5 A Corner Body 
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Man»« !R Coalmenm of Flexible end Body ll.irmnnlri 
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